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1. Introduction 


The apparent matter-antimatter asymmetry of the universe remains an outstanding mystery of 
physics that cannot be explained within the standard model (SM) [jj. This disagreement between 
the SM and observation motivates the existence of beyond the standard model (BSM) theories with 
additional violation of CP and baryon number {B) conservation. Many BSM models with additional 
symmetry violation can be constructed, and experimental guidance is necessary to constrain the 
space of viable BSM models and determine what theory describes our matter-dominated universe. 

Experimental measurements of, or bounds on, baryon number violation rates can only be di- 
rectly compared to theoretical predictions of transition rates between hadronic States. Reliable 
calculations of hadronic observables must include the strong nuclear forces present, and perform- 
ing top-down calculations of hadronic transition rates for every BSM model of interest is highly 
impractical. Instead, baryon number violating processes can be treated in a model-independent 
way by considering the standard model as an effective field theory that includes baryon number 
violating operators composed only of quarks and other SM degrees of freedom. The predictions 
of particular BSM models of interest can be straightforwardly parametrized in terms of local SM 
operators. Once the model-independent matrix elements of these operators have been determined, 
BSM models can be directly confronted with experimental constraints. 

The most relevant baryon number violating operator that can be added to the standard model 
is a dimension 6 operator that permits proton decay. Experimental constraints on proton decay 
rates place stringent bounds on matrix elements of this operator and therefore the scale of NB = 1 
physics. Many phenomenologically attractive BSM theories include a suppression of these well- 
constrained NB = 1 effects, see Ref. 0 for a review. In these theories, the baryon number violating 
effects most accessible in low-energy experiments are typically NB = 2 neutron-antineutron transi- 
tions. Neutron decay searches have a long history [^, hut extracting neutron transition rates from 
observations of nuclei remains challenging. Experiments using reactor neutron beams can pro¬ 
vide clean measurements of the n — n transition rate 0, and current experimental bounds have the 
potential to be greatly improved with new measurements |^. 

Tuming precise experimental data on n — n transition rates into precise constraints on BSM 
physics presents theoretical challenges. Naive dimensional analysis alone cannot be used to reliably 
estimate transition rates-varying the value of Aqcd used for naive dimensional analysis of a six- 
quark operator matrix element leads to a thousand-fold transition rate variation. The requisite 
nuclear matrix elements have also been estimated within the MIT bag model 0, hut there is no 
way to determine the theoretical uncertainty present. 

The only technique currently available for turning experimental n — n transition rate data 
into constraints on BSM physics with quantified uncertainties is lattice quantum chromodynam- 
ics (EQCD). Preliminary EQCD calculations of the required matrix elements have been performed 
1^]. These simulations include large systemic uncertainties arising from the need to reläte lattice 
regularized matrix elements defined at GeV scales to MS renormalized matrix elements defined 
at high scales where reliable BSM predictions have been made. This report describes an ongoing 
calculation of these operator renormalization effects. 
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2. Operator Renormalization 


The Hamiltonian for « — TT transitions can be expanded in a complete basis of six-quark opera¬ 
tors, explicitly constructed below. Denoting these operators by ^ the Wilson coefficients predicted 
by a particular BSM theory by C“, the Hamiltonian density for n — TT transitions is expanded as 




( 2 . 1 ) 


where the first equality involves bare operators and coefficients and the second renormalized oper¬ 
ators and coefficients defined at an arbitrary renormalization scale jJ.. These are related by renor¬ 
malization coefficients that are defined for a particular renormalization scheme and scale as 
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(M)=^ 


ab 

scheme 


( 2 . 2 ) 


BSM models can typically be used to perturbatively calculate Wilson coefficients of MS renormal¬ 
ized operators at high scales. Model independent relations between these coefficients and measur- 
able n — n transition rates can only be provided once lattice regularized matrix elements extracted 
from LQCD simulations are related to MS renormalized matrix elements at high scales. 

Since the MS renormalization prescription can only be applied to dimensionally regularized 
(dim reg) matrix elements, relating lattice and MS matrix elements necessitates the introduction 
of an intermediate renormalization scheme that can be applied to both lattice and dim reg matrix 
elements. The Regularization Invariant Momentum (RI-MOM) scheme is well suited for this |||], 
and has been utilized for practical operator renormalization of kaon oscillation ||^ and proton decay 


OlOO matrix elements. 

The RI-MOM renormalization scheme, explicitly defined below, is implemented by demand- 
ing that vertex functions calculated at a reference momentum po are fixed to a particular value. This 
definition can be applied non-perturbatively to LQCD matrix elements in order to determine the 
coefficient matrix Z]j““(po) relating lattice regularized and RI-MOM renormalized operators. The 
RI-MOM scheme can also be applied perturbatively to dimensionally regularized matrix elements 
to determine Z^^”^(po) and necessarily removes all poles in D — 4 to produce fini te renormalized 
matrix elements. This means that the ratio of Z^/“{pq) to Zj^{pq) is a finite quantity that can be 
reliably calculated in perturbation theory provided »^(po) "C 1 . 

Combining the matching factor given by this ratio with standard renormalization group analy¬ 
sis, the relation between MS operators at arbitrary scales and lattice operators with quark momenta 
Po is given by 
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(2.3) 


where jSo and jSi are the one-loop and two-loop coefficients of the QCD p -function and 70 and 71 
are the one-loop and two-loop coefficients of the AB = 2 operator anomalous dimension matrix. 

Keeping both discretization and perturbative errors parametrically suppressed requires the hi- 
erarchy of scales Aqcd <C po “C where a is the LQCD lattice spacing. BSM matching calcu- 
lations are performed at much higher scales than present inverse lattice spacings, and so Eq. 2.3 
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will be generally applied to cases where <C ^^(po)- In this case both one-loop contributions 

to the matching factor ( ) nnd the two-loop anomalous dimension 71 provide 0{as{po)) 

operator renormalization effects. The one-loop anomalous dimension has been known for some 
time [ 0 > but both the one-loop matching factor and two-loop anomalous dimensions are unknown 
and form the targets of our investigation. 


3. Operators and Vertex Functions 


Six-quark operators that are Lorentz invariant, local, standard model gauge invariant, and have 
the correct flavor structure to describe n — n transitions can be expressed in a basis of operators of 
the form' 

^X\X2Xi j ^^X2^l)^^mP'^X2^ri)T{ij}{kl}{mn}i (3-1) 

^X\X2X2 

^XiX2X3 ~ ^^X\^j)^^k^^X2^l)i.^m^^X3^r)'^[ij\[kt\{mn}i 

where •'' = are chirality labels with = 3(1 ± 75 ), C is the charge conjugation matrix, 

and above are color labels. The only tensors that can combine six quarks of this flavor 

structure into color singlets are 


'^{ij}{kl}{mn} — ^ikm^jln ^jkm^iln ^ilm^jkn ^ikn^jlm: (3-2) 

'^[ij][kl]{mn} ~ ^ijm^kln T ^ijn^klmi 

where Sijk is the Levi-Civita symbol. These tensors are symmetric {/ j} or antisymmetric [/ j] in each 
pair of indices as shown. Operators with more complicated Dirac matrix structures are not indepen- 
dent and can be expressed as linear combinations of those above through Fierz transformations. The 
trivial constraints and non-trivial constraints ^ ^^xxx' 

immediately reduce the number of independent operators of this form from 24 to 14. Parity further 
reduces the number of independent operators to 7. Chiral symmetry provides further constraints, 
but maintaining a complete operator basis when utilizing chiral symmetry breaking lattice regular- 
izations requires a basis of 7 independent operators such as 

^LRL-: ^RRR-: ^LLRi ^LRR^ ^RRr}- (3-3) 


Vertex functions for each operator can be constructed (in both perturbation theory and LQCD) 
by Wick contracting with neutron and antineutron interpolating operators 


(/.a\C(PYSria 
\ njklmn 




ufd[di 


amp 


(3.4) 


where the subscript amp refers to the prescription of amputating extemal legs by taking 




— l\a'a 
q )i'i 


(3.5) 


*We thank Brian Tiburzi and Sergey Syritsyn for helpful discussions on the symmetry transformation properties of 
six-quark operators and optimal operator basis construction. 
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where Sq is the propagator for quark flavor q. The tree-level vertex functions for these operators 
are given explicitly by 

^ „ ^ p'Y,s',r]'W\^a'l3'YS'r]'a'rj. .o 

■^(0 )) ■ j^i ^jid^idp d^,d^i ^2 "^3 d^{i'j'}{k'l'}{m'n'} (3-6) 

-^^{ij}{kl}{mn}[^\ ^2 ^3 +^1 ^2 ^3 \ {ij}{km}{ln} [^l ^2 ^3 +^1 ^2 ^3 J 

+^T{,j}^kn}{lm} +XfX^^X^^] , 

^"o))lZr = “ (^dld^Ji'dlu]u^t^\^^ (3.7) 

- ^-l(ik)(jl){mn}[Xx X^ A 3 +ÄJ Ä 2 A 3 J + 27[Aj A 2 A 3 +Aj A 2 A 3 J 
+^^(ik)(jn){lm}\XT^ A 2 A 3 +Aj A 2 A 3 \ + Il qiYjk){mn}[X\ X 2 A 3 +Aj A 2 A 3 J 

+^l(il){jm){kn}[^\ ^2 "^3 ^2 ^3 \ + ^1 (il){jn){km}[^\ ^2 "^3 ^2 ^3 J 

+^l{im){jk){In}[Xi ^2 A 3 +Aj A 2 A 3 J + 27 (/„j)( 27 ){^„} [Aj A 2 A 3 +Aj A 2 A 3 J 

_1_9T fy^^y^^^y -1-y^^y®*^ y^^l _L 9T [yCCOyPYySTI , v/t/yggv^m 

+^l{im){jn){kl}[Xi X 2 A 3 +Aj A 2 A 3 J + 27 (,-„)( 27 .)p,„} [Aj A 2 A 3 +Aj A 2 A 3 J 




_i_9T Ty^‘^yP®y-1-yP° y'^‘^y_L 9T [yaaypiiyYQ i yP^yaayy 

where X,- = CP^^, {ij) denotes {/ j} for and [/ j] for ^A, and we have suppressed chiral indices on 
‘‘A. 

The RI-MOM renormalization condition is formally defined by a condition on these vertex 
functions evaluated with all quarks at a reference (Euclidean) momentum po and 

-3 r^h^aPYSna , 

' ijklmn 




aPYSrja 
i jklmn 




= 5“* 


(3.8) 


where Z‘^qj^ is the RI-MOM wavefunction renormalization factor defined in Ref [^, is 

the matrix of renormalization coefficients appearing in Eqs. (^), (^), and projects onto the 
tree-level vertex function A“qj according to the defining relation 

,aPYSna(.x, 


{py. 


i jklmn 


A' 


(0)7 


i jklmn 


= 5“”. 


(3.9) 


There is some freedom in constructing projectors that satisfy Eq. (^^. The set we have constructed 
is given explicitly by 

apYSa 1 




{Plrl) , 


i jklmn 

aPySa 
i jklmn 

aPyScj 
i jklmn 


55296 
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{U}W}{>nn}^R 

r C(S yYP yV^ 


(3.10) 
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Projectors satisfying Eq. (3.S) are useful for decomposing matrix elements into a desired oper¬ 
ator basis in more general contexts than enforcing the RI-MOM condition Eq. (^). In particular, 
we can express the MS renormalization condition as 


7^(7‘l \ ^ ^ ^aj5rSr]a ^ ^ac 

^MS\^lÄs) Jijklmn \ PoleJijkbnn ^ ’ 


(3.11) 


where denotes the dim reg vertex function contributions proportional to l/e = 2/(4 — D) — 
7 £ -|-ln 47 r, Z^is the usual MS quark wavefunction renormalization factor, and we have suppressed 
dependence on the renormalization scale jj.. 

Using projectors to decompose matrix elements into a particular operator basis is essential for 
EQCD relations where the analytic relations between full and tree-level vertex functions are un- 
known. Given the complicated tensor structures of our vertex functions, the projector definitions of 
Zri and have practical advantages for our perturbative calculations as well. Most importantly, 
applying the projector definitions greatly simplifies the tensor decomposition of the many spin, 
color, and flavor structures produced by one-loop and especially two-loop diagrams in dim reg. 


4. Perturbative Calculation 


A calculation of the one-loop RI-MOM/MS matching factor includes 15 one-loop Eeynman 
diagrams. These can be straightforwardly expressed as contractions of four independent tensor 
integrals, two for Eeynman gauge and two more for general gauge needed for consistency with 
EQCD simulations performed in Eandau gauge. 

The contractions of these tensor integrals with the Dirac structures of the diagrams introduces 


subtleties, however. Applying the projectors to structures such as i 
introduces traces of 4 7 ^ times 75 that are ill-defined in naive dimensional regularization. While 
these structures can be related to our original operator basis through Eierz transformations, these 
Eierz transforms are only valid in four dimensions. To consistently account for the fact that our 
operator basis is only complete in D = 4 we follow the approach of Ref. 10 and introduce evanes- 
cent operators proportional to D — 4 that complete our basis in D dimensions. Both the one-loop 
matching factors and two-loop anomalous dimensions depend on the basis of evanescent operators 
chosen and are ambiguous until this basis is specified. 

320 two-loop diagrams contribute to the MS anomalous dimension matrix. Standard multiloop 
techniques can be used to perform tensor reductions of the independent tensor integrals appearing 
and to use integration by parts algorithms to recursively decompose the resulting scalar integrals 
into a minimal set that can be evaluated straightforwardly [|T^, M, 15]. Slightly more than half 
of the independent two-loop topologies can be treated as simple generalizations of the diagram 
topologies appearing in four-quark operator renormalization calculations, but the remaining dia¬ 
grams involving gluon exchange between all three spin-singlet diquark structures in our operators 
have no direct four-quark analogs. Proper treatment of evanescent operators in these diagrams 
in particular requires the introduction of new evanescent operator structures that do not appear in 
four-quark operator calculations. 




rs t}T]o 
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5. Conclusion 

Extracting predictions of « — «transition rates with quantified uncertainties from BSM theories 
relevant for baryogensis requires LQCD simulation as well as perturbative operator renormalization 
analysis. Aeeurately ealeulating relations between lattiee regularized and MS renormalized matrix 
elements requires both non-perturbative renormalization of LQCD matrix elements in a regular- 
ization invariant seheme sueh as RI-MOM and perturbative matehing between RI-MOM and MS 
renormalization faetors. In addition to this one-loop renormalization seheme matehing ealeula- 
tion, the two-loop MS anomalous dimension matrix is needed for aeeurate relations between BSM 
matehing ealeulations and LQCD simulations at different seales. 

We have begun an ongoing ealeulation of both of these perturbative quantities using suit- 
ably generalized four-quark operator renormalization teehniques as well as operator projeetors that 
simplify two-loop tensor deeompositions. Completing this ealeulation will remove a signifieant 
systematie uneertainty in using n — n transition rate measurements to eonstrain BSM physies. Our 
two-loop ealeulation is well underway, and we will publish results on both the one-loop matehing 
faetors and two-loop anomalous dimension matrix onee the two-loop ealeulation is eomplete. 
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